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ABSTRACT 

We present high-resolution radio continuum observations of the nascent starburst in the metal- 
poor galaxy SBS 0335— 052. These radio data were taken with the Very Large Array and include 
observations at 0.7 cm, 1.3 cm, 2 cm, 3.6 cm, and 6 cm. These observations enable us to probe 
the thermal radio nebulae associated with the extremely young star-forming regions in this galaxy. 
Two discrete and luminous star-forming regions are detected in the south of the galaxy that appear 
to be associated with massive star clusters previously identified at optical wavelengths. However, 
the remaining optically-identified massive star clusters are not clearly associated with radio emission 
(either thermal or non-thermal) down to the sensitivity limits of these radio data. The spectral energy 
distributions of the two radio-detected clusters are consistent with being purely thermal, and the entire 
region has an inferred ionizing flux of ~ 1.2 X 10^3s~\ which is equivalent to ~ 12,000 "typical" O- 
type stars (type 07.5 V). The observations presented here have resolved out a significant contribution 
from diffuse non-thermal emission detected previously, implying a previous episode of significant star 
formation. The current star formation rate (SFR) for this southern region alone is ^ 1.3 MQyr~^, 
or ~ 23 M0yr~^kpc~^, which is nearing the maximum starburst intensity limit. This SFR derived 
from thermal radio emission also suggests that previous optical recombination line studies are not 
detecting a significant fraction of the current star formation in SBS 0335— 052. From model fits to 
the radio spectral energy distribution, we infer a global mean density in the two youngest clusters 
of Ue > 10"^ — lO'' cm~'^. In addition, a comparison between the compact and diffuse radio emission 
indicates that up to ^ 50% of the ionizing fiux could be leaking out of the compact Hii regions; this 
in is agreement with previous work which suggests that the interstellar medium surrounding the natal 
clusters in SBS 0335—052 is porous and clumpy. 

Subject headings: galaxies:individual(SBS 0335-052) - galaxies: star clusters - galaxies: starburst — 
HII regions - stars: formation 
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1. INTRODUCTION 

Super star clusters (SSCs) are the most massive and 
dense type of young stellar clusters, and represent the 
most extreme mode of star formation in the local uni- 
verse. S SCs are common in vigorous star f ormation 
episodes (jHunter et al.l 119941 : iKeto et al.ll2005f ). appar- 
ently requiring high pressure and interstellar gas column 
densi ty for their formation (e.g., Elmcgrcen & Efremo3 
ll997t[Biilett et al.ll2002f ). A large amount of evidence col- 
lected since the launch of the Hubble Space Telescope sug- 
gests that at least some SSCs are the progenitors of glob- 
ular clusters; they have similar masses, sizes, and stellar 
density, and evolutionary models suggest that evolved 
SSCs, if gravitationally bound, woul d have properties 
very similar to pre s ent-dav globulars (iBarth et al.|[l995t 
iMeurer et al]ll995t IHo fc Finppenkolll996l : iLarsen et all 
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I200H IWhitoorillOOl). However, despite these strong 
links between globular clusters and SSCs, we do not know 
how the extremely low metallicity in the early universe 
affected massive star cluster formation. A low metal con- 
tent may affect star formation in a variety of ways, in- 
cluding cooling and pressure in the birth cloud, hardness 
of the stella r spectra, and e volution of the stars them- 
selves (a^gJSchaere3j200^JSm^ & CrowtheJ 
I2002t iTumhnson. Venkatesan. fc Shullli2004i : iBatQi2005f) . 
One way to approach this question is to study SSC for- 
mation in extremely low-metallicity systems in the local 
universe. 

SBS 0335—052 is a unique blue compact dwarf galaxy 
(BCD) in the nearby universe because of its very high 
star formation rate (SFR) and its extremely low metal- 
licity; it is thus ideal for our study. It was discov- 
ered i n the Second Byurakan Survey by llzotov et al.l 
(|1990( ) who noted its particularly low oxygen abun- 
dance of ~ 1 2 -I- log(0/H) = 7.3, confirmed by subse- 
quent studies dMelnick. Hevdari-Malaveri. fc Leisvl[l992l : 
llzotov et ahl 1200 It ). SBS 0335-052 is the most metal- 
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Fig. 1. — Multicolor Hubble Space Telescope ACS image of SBS 0335-052 constructed us ing the Ha FR656N filter (red), the visual 
F550M fiher (green), and the uhraviolet F220W+F330W filters (blue). The SSCs identified bv lThuan eTall (IT997I ) are labeled. The image 
is shown in the standard orientation with North being the top of the image and is ~ 10" X 10". 

iPapaderos et al.lll998^ . Toward the north, the filamen- 
tary and irregular structure of the emission suggests a 
shell carved out by supernovae (SNe). The SSC s (1- 
6, according to the notation of iThuan et al.l [19971 ) are 
roughly aligned in the southeast to northwest direction. 
iReines. Johnson, fc HuntI ()2008f ) confirm an age gradi- 
ent of the clusters across the galaxy (with the northern- 
most clusters being the oldest) that is consistent with a 
large-scale disturbance crossing the galaxy at a speed of 
~ 35 km s~^ and triggering star formation en route. 

Although the optical colors of the two youngest 
SSCs 1-1-2 are relat ively blue, th ey are also strong in- 
frared (IR) emitt e rs (iThuan et al, 1999: Hunt et al. 200j; 
Dale et all [Mil iPlante fc Sauvagd l2002t iHouck et alj 
2004 IReines. Johnson, fc HuntI \20m . Extinction es- 
timates around these two clusters vary greatly de- 
pending on the wavelengths of the observations 
used. Low extinctions (Ay ^ 0.5 mag) arc de- 
rived from optical observations (e.g.. Ilzotov et al.lll997l : 
IReines. Johnson, fc HuntI 120081 ) whereas mid-IR obser- 
vations indicate high extinction (Ay > 12 mag) regions 
surrounding the young clusters (e.g.. lThuan et anil999l : 



poor galaxy known with a SFR ^1 Mp, yr ^, > 20 times 
higher than in I Zw 18 (jFanclli. O'ConnelL fc ThuanI 
[l988t iHunt eTall l2005bD . SBS 0335-052 is associated 
with an Hi condensation within the same large Hi 
cloud (64kpcx24kpc) as SBS 0335-052W, which is lo- 
cated at a second H i peak about ~22kpc to the west 
(jPustilnik et al.ll200l f). The neutral gas in SBS 0335—052 
is evidently not pristine, having ro ughly the same oxy - 
gen abundance as the ionized gas (jThuan et al.ll2005r ). 
Analyses of the broadband colors and spectral models 
of SBS 0335— 052 suggest that it may be a local ana- 
log to primordial star forming sites in the early u niverse 
Jriiuan et al. 1997; Izotov ct al. 1997; Papader os et al.l 
[l998HPustilnik et al.ll2004D . 

The star formation in SBS 0335—052 is concentrated 
mainly in six SSCs^ (shown in Fig. [H), surrounded 
by a very blue underlying envelope f|Thuan et al.lll997l : 

^ Based on observations made with the Hubble Space Telescope, 
obtained from the Data Archive at the Space Telescope. Science 
Institute, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under NASA contract NAS 5-26555. 
These observations are associated with program #10575. 
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Plante fc Sauvagell2002t iHunt et al ll2005at iHouck et all 
20041 ). iReines. Johnson, fc HuntI (|2008l ) find evidence 
for a porous and clumpy ISM surrounding the young 
clusters, which can naturally account for the apparently 
discrepant extinction estimates found in the literature. 
Dense dust clumps heated by the impinging UV stellar 
continuum provide the high extinction estimates from 
mid-IR observations, whereas diffuse inter-clump regions 
sampled by optical observations yield low measured ex- 
tinctions. 

Previous low-resolution radio observations of 
SBS 0335— 052 showed free-free absorption on a 
global scale , and a significant non-thermal component 
([Hunt et al.l 120 0^ . However, the beam size of those 
images was insufficient to disentangle the spatial distri- 
bution of the thermal and non-thermal emission. In this 
paper, the radio continuum properties of SBS 0335—052 
are re-examined at significantly higher resolution with 
the goal of probing individual star-forming regions. We 
present new high-resolution radio continuum observa- 
tions of SBS 0335—052 to better separate the different 
emission mechanisms spatially, and probe the thermal 
radio nebulae. We describe our observations in in §2, 
discuss the morphology of the natal clusters in §3, and 
analyze the physical properties of the clusters in §4. The 
relationship between the clusters and the overall galaxy 
are presented in §5, and §6 provides a general summary 
of the results. We adopt a distance of 55.7 Mpc to 
SBS 0335-052, which assumes i/o = TOkms"! Mpc'^, 
an d correction to the CMB reference frame as described 
in iHunt et"aLl ()2005af ). This corresponds to a spatial 
scale of 270 pc arcsec"^. 



High-resolution radio observations of SBS 0335— 052 
were obtained with the Very Large Array (VLA)^ from 
2003 June to 2005 August. Observations were obtained 
at C-band (5 GHz, 6 cm), X-band (8 GHz, 3.6 cm), U- 
band (15 GHz, 2 cm), K-band (22 GHz, 1.3 cm), and 
Q-band (43 GHz, 0.7 cm); these are summarized in Ta- 
ble [1] The high-frequency observations at K-band and 
Q-band utilized fast-switching to a nearby phase calibra- 
tor with cycle times of ~ 2 minutes in order to mitigate 
the effect of atmospheric changes. 

These radio data were reduced and calibrated using the 
Astronomical Image Processing System (AIPS). When 
available, models of the flux calibrators were used in or- 
der to exploit the full uv coverage of these sources. Most 
of the data sets included observations of more than one 
flux calibrator in order to cross check the calibration solu- 
tions. We estimate that the absolute uncertainty in flux 
calibration is < 5%, based on scatter in the VLA Flux 
Density Calibrator database for the calibrators used in 
this program. 

Because our aims depend on relative flux densities at 
different frequencies, care was taken to obtain the best- 
matched uv coverage possible at each frequency. Each 
frequency was observed in two separate array configu- 
rations chosen to obtain relatively well-matched synthe- 
sized beams. Additionally, all of the combined data sets 
were trimmed to have an identical minimum uv spacing 
of 30 kA in order to roughly match their sensitivity to ex- 
tended structure, and the "robust" parameter was varied 
between -1 and 1 (slightly uniform to slightly natural) in 
order to obtain better matched synthesized beams. Nev- 
ertheless, the uv coverage at each frequency cannot be 



2. OBSERVATIONS 



The National Radio Astronomy Observatory is a facility of 
the National Science Foundation operated under cooperative agree- 
ment by Associated Universities, Inc. 
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TABLE 2 

Imaging Parameters 
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TABLE 3 
Radio Flux Densities 





^1.3cm 


-Fb.Ocm 


^3. 6cm 


-Ffi.Ocm 


Source 


(xlO-4 Jy) 


(xlO-4 Jy) 


(xlO-4 Jy) 
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81 


1.7 ±0.3 


1.8 ±0.5 


1.8 ±0.4 


1.4 ±0.2 


S2 


1.2 ±0.2 


1.0 ±0.2 


0.8 ±0.2 


0.8 ±0.2 


Region S"* 


5.6 ±0.6 


5.8 ±0.6 


4.6 ±0.5 


3.4 ±0.4 



^ Region S corresponds to an aperture of 0.5 radius centered on tlie radio 
emitting region that encompasses SI and S2. 



perfectly matched, and thus there may be shght varia- 
tions between the sensitivity to different spatial scales 
at each frequency. An additional set of images was also 
made using the greatest possible sensitivity at each fre- 
quency using purely natural weighting (robust=5) and 
no restrictions in uv coverage. The resulting imaging 
parameters are listed in Table [2l 

Unfortunately, the highest frequency observations at 
Q-band were consistently plagued by bad weather and 
instrumental problems. As a result the rms noise in the 
combined Q-band data set is too high to be useful for 
the purposes of this project, and is not analyzed further 
in this paper. However, in the final data sets obtained 
in August 2005, a luminous transient object appeared, 
which will be discussed in a separate paper. 

Determining the flux densities of non-point-like sources 
observed with an interferometer amid complex back- 
grounds is notoriously prone to large uncertainties. In 
order to obtain the most accurate values possible, we 
employed several techniques to measure the flirx den- 
sities. First, the peak flux densities were measured in 
mJy/beam; if the objects are point-like, this is equal to 
their total flux density in mJy. Second, aperture pho- 
tometry was performed on each wavelength using a va- 
riety of identical (irregular) apertures and annuli using 
the AIPS++ viewer tools. Finally, the sources were fit 
with Gaussian profiles using the AIPS task JMFIT. The 
quoted flux densities and uncertainties in Table [3] reflect 
these results. In addition, photometry was performed on 
the high-sensitivity images using a 0.5" radius aperture 
(= 135 pc) centered on the entire region of radio emission 
"region S" , and these results are also quoted in Table [31 

3. LOCATION AND MORPHOLOGY OF THE RADIO 
EMISSION 

The compact radio emission in SBS 0335— 052 is re- 
solved into two dominant sources SI and S2 located in 
the southern region of the galaxy that appear to be asso- 



ciated with two SSCs 1 and 2 identified bv lThuan et al.l 
(|1997f ). The radio emission at 1.3 cm, 2 cm, 3.6 cm, and 
6 cm is shown overlaid on an HST/ AGS F550M image in 
Figure [2l 

The exact location of the radio sources relative to the 
optical SSCs is difficult to pinpoint precisely because of 
the uncertainty in the HST astrometry. The astrometry 
for the VLA images is accurate to within ^ O'.'l and is 
considered absolute by co mparison. We starte d first with 
the same HST image as in iHunt et al.l (|200l . astromet- 
rically calibrated with five stars in the USNOA2.0. The 
astrometric rms uncertainty is 0'.'54, roughly twice the 
size of the highest resolution synthesized beams of the 
radio data presented here. Using the nominal astrome- 
try, the southernmost radio peak lies ~0'.'4 south of the 
two brightest SSCs. 

We then attempted to improve the relative pointing 
uncertainty between the optical and radio images by 
bootstrappi ng from the morphology of the P act emission 
presented in lReines. Johnson, fc Hun^ ()2008f ) . Using this 
method, we first matched the continuum-subtracted Paa 
emission to the radio, and applied that astrometric solu- 
tion to the non-continuum-subtracted Paa image, which 
has many features in common with the optical image. 
Subsequently, the optical image was aligned with the 
non-continuum subtracted Paa image to obtain the fi- 
nal astrometry. The shift in the astrometry relative to 
the USNO solution is '^0'.'4, within the expected uncer- 
tainty of the absolute positions. The resulting alignment 
of the optical image clearly associates the radio sources 
SI and S2 with the optically identified SSCs 1 and 2 as 
shown in Figure [2l Given this apparent correspondence 
between the optical and radio sources, hereafter we refer 
to SI and S2 as SSCs 1 and 2. 

None of the other SSCs identified optically by Thuan 
et al. have firmly-detected radio emission, although 
SSCs 5 and 6 have possible marginal detections in the 
1.3 cm K-band image. In addition, there is a string of 
K-band ^ 3 — 5cr peaks extended toward the north of 
SBS 0335-052 (see Fig. [2|). Individually, any one of 
these peaks would not necessarily be considered a source; 
however, the spatial correlation is suggestive of a more 
extended feature that is largely resolved out by these 
high-resolution observations. We hypothesize that this 
fea ture may be associa ted with the supershell observed 
bv iThuan et al.l ()1997f ) and/or an outflow extending in 
that direction. 

The flux density at high frequencies observed for Re- 
gion S in this p aper compared to the previous work at 
lower resolution (i Hunt et al.ll2004 ) suggests that even in 
these high-resolution observations the spatial filtering in- 
herent in array observations is not "resolving out" any of 
the thermal flux. At 1.3 cm and 2 cm, our high-resolution 
observatio ns recover all the flux in SBS 0335— 052 re- 
ported bv lHunt et al.l (|2004l ). This implies that the size 
of the entire thermal source at high frequencies must be 
smaller than or comparable to the 0.5"aperture used to 
measure the flux density of Region S. 

In contrast, the 2 cm radio emissio n at a resolution 
of I'.'e presented by iHunt et al.l (|2004f ) is peaked -0'.'5 
(~100pc) north of sources SSCs 1 and 2. Fig. [3] illus- 
trates the overlay of the 20 cm emission on the 1.3 cm 
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Fig. 2. — VLA contours overlaid on an HST ACS F550M image in gray-scale, (a) : VLA 6 cm contour levels of 3, 4, 5, 6 X cr (18 ^Jy/beam). 
(b): VLA 3.6 cm contour levels of 3, 4, 5, 6, 7 X o- (12 /^Jy/beam). (c): VLA 2 cm contour levels of 3, 4, 5, 6 X cr (17 /^Jy/beam). (d): VLA 
1.3 cm contour levels of 3, 4, 5, 7, 9 X cr (11 ^Jy/beam). 



map^. Because the low- frequency radio flux is domi- 
nated by synchrotron emission, this non-thermal emis- 
sion is likely to be associated with SNe from a distinct 
episode of star formation with an age of >3.5Myr. At 
6 cm, the lowest frequency of the observations presented 
he re, we recov e r only 44% of the global flux measured 
by iHunt et~an (|2004[ ). This would imply that the non- 
thermal emission is much more diffuse than the ther- 
mal; indeed if the s ignificant synchrotron halo of I Zw 18 
()Hunt et al.l l2005"bl ) . another low-metallicity BCD, were 
placed at the distance of SBS 0335— 052, it would be re- 

T he positional discrepancy between this work and IHunt et al.l 
II2004I ) is a consequence of the new HST images and astrometry 
presented in this paper. 



solved out by our observations (see H5A\\ . 

4. PHYSICAL PROPERTIES OF THE NATAL BSCS 

4.1. Modeling the Spectral Energy Distributions 

To gain insight into the physical conditions in the ra- 
dio sources present in SBS 0335— 052 we apply a general 
ther mal-|-non-the r mal fi t with an absorption component 
as in IHunt et al.1 (|2004D . These are simple models that 
assume a uniform-density homogeneous ionized medium. 
The realistic physical conditions in the Hii regions are 
likely to be more complex, including density variations 
and deviations from a sphere, but the available data do 
not support more complicated models. We can express 
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Fig. 3.— L-band (20 cm) contours from lHunt et al.l 1I2OOI ') overlaid on (left) the K-band (1.3 cm) image of SBS 0335-052 and (right) the 
ACS F550M image of SBS 0335-052. Although at lower resolution, the centroid of the 20 cm source is offset from the thermal emission by 
~ 0'.'5 to the north. 



the thermal radio free-free absorption coefficient 



~ 0.08235 ( ^ 



-1.35 



GHz 



-2.1 



pc 



where rie is the electron density of the ionized gas, T is 
the ionized gas temperature, and v is the frequency. This 
means that the free- free optical depth Tg — L (i^path 
length through the region) is: 



Tff ~ 0.08235 



-1.35 



EM 



pc cm 



GHz 



-2.1 



(2) 



where EM is the linear emission measure ('^ 

By virtue of Kirchhoff's law for thermal emission, 
which relates and the radio volume emissivity , we 
write the thermal radio flux as a function of EM: 



C--= ^(5.95x10-^) 



f_EM^)(^)2 mJy 

y pc cm / V arcsec / 



(3) 



(j) is dimensionless and depends on the geometry, equal 
to 7r/6 for a spherical region of constant density and di- 
ameter (FWHM) 9, and it / A for a cylindrical region o f 
diameter and length 9 (e.g., Mczgcr & Henderson 1967). 
(j) also accommodates a filling factor, which here is as- 
sumed equal to unity. We take into account the optical 
thickness of the radio emission through assuming a 
geometry in which the emission and absorption are in- 
termixed in a homogeneous medium. The thermal radio 
flux can then be expressed as: 



J V 



l-CXp(-Tff) 



thil 



= [1 - exp(-Tff)] (7.225 x IQ-^) 
X i^-)^ mJy 



(k) (gHz) 



(4) 



Equation 0] shows that in the limit of optically-thin ra- 
dio emission where the absorption term [1 — exp(— Tff )] oc 



th 



-0.1 



Alternatively, when Tff 3> 1, the ra- 



dio emission is completely optically thick, and / 



Hence we are not making any assumptions a priori about 
the nature of the thermal radio emission^. A spheri- 
cal geometry is assumed for the thermal emission, and 
the electron temp erature Tg was taken to be 20000 K 
(jlzotov et al.lll999D . 

Allowing for the most general models that also in- 
clude non-thermal emission introduces two additional 
free parameters (the non-thermal flux density and spec- 
tral slope), leaving the models unconstrained. Neverthe- 
less, for completeness we also modeled these regions with 
the inclusion of a non-thermal component to assess the 
degree to which the results might be affected. We assume 
a screen geometry for the free-free absorption and write 
for the non-thermal radio flux /"*: 



= exp(-Tff) n 



mJy 



(5) 



where /"^' is the non-thermal (unabsorbed) flux at fre- 
quency vq. The total flux in this case at given frequency 
V is the sum of Equations [3] and [51 

The data were fit to models allowing for both thermal 
and non-thermal contributions using a minimization 
algorithm, and contour values of the reduced values 
and resulting spectral ener gy distributions ( SEDs) are 
shown in Figures[4^-c. As in lHunt et al.l (|2004l ). we fit the 
data using three free parameters: /"*, diameter 0, and 
emission measure EM. We ran models separately for val- 
ues of ant ranging from —0.45 to —0.8, which would en- 
able us to judge, albeit crudely, the best-fit non-thermal 
index from the lowest x^- The best- fit parameters de- 
rived from these models are listed in Table S] and dis- 
cussed in the following. In each case, the radio fiux den- 
sities are consistent with a purely thermal origin. 

4.2. Radii and Densities 

The best-fit models with only thermal emission are 
shown as solid curves in Fig. [4l The models also allow 
for an additional non-thermal contribution, although this 
only resulted in a marginally better fit for for SSC 1 and 
is shown as a dashed line. The best-fit parameters and 



These are the same models as in IHunt et all 1I2OOI 1. but ex- 
pressed in a more general way. 
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TABLE 4 

Properties of the Thermal Radio Sources Inferred from the 
Spectral Fits 









lie 


Q Lyc 




Source 




(pc) 


(10^ cm-^) 


(1052 s-l) 


(10^ Mq) d 


SI 


0.2 


5.0 


4.3 


3.5 


0.5 


Slnt ° 


0.1 


4.5 


4.8 


3.5 


0.5 


S2 


0.6 


2.8 


7.9 


2.5 


0.4 


Region S 


0.3 


7.9 


3.2 


11.7 


1.7 



It is clear from the reduced values that these fits are not well- 
constrained. See § 14.21 for discussion.'' See Figure |4] for range of 
acceptable values.'^ Qhya values as determined from the 1.3 cm flux 
densities, assuming an Hll region temperature of 20 000 K. Adopt- 
ing t he models of SB99 for a burst at age IMyr IILeitherer et al.l 
119991) .° Source SSC 1 had a marginally better fit for models that 
included a small amount of non-thermal emission. 

the physical properties inferred from these models are 
reported in Table H) However, as can be seen in Fig- 
ures H^-c, it must be stressed that in all cases presented 
here, there is a range of parameters that fit the data 
with nearly the same goodness-of-fit. It is important to 
understand the degeneracy inherent in fitting thermal 
radio SEDs of unresolved sources. Specifically, at high 
frequencies where the resolution is sufficient to separate 
star-forming regions, sources also tend to be optically 
thin or only slightly optically-thick. In this regime, ra- 
dius and density are essentially degenerate. From the 
plots of the reduced minima shown in Figures U^-c, it 
is clear that a range of values for radius and density are 
consistent with the data in hand; thus the "best-fit" val- 
ues quoted in Table 0] should not be strictly interpreted. 

The radii of the two thermal sources associated with 
SSCs 1 and 2 inferred from our fits are ~ 3 — 6 pc (see 
Table |3]). These sizes are sev eral times smaller than 
the radii of 8-15 pc deduced bv iHunt et alj (|2004[) from 
lower-resolution observations, but typical of SSCs. They 
are also considerably smaller than those measured by 
[Thompson et al. (2008) with NICMOS at 1.9 /im (16- 
18 pc). However, the diffraction limit of HST at this 
wavelength of ~0.2arcsec (54 pc) makes direct measure- 
ments with NICMOS potentially unreliable. Likewise, 
the radii inferred here are not directly measured, and 
therefore these sizes should be considered estimates. 
Moreover, given the youth of these clusters, they are 
unlikely to be dynamically relaxed into simple spheri- 
cal systems, but rather complex structures that are not 
resolved by any of the available data. 

The inferred EMs of 2 — 4 x 10® pc cm~^ are about 
10 times higher than that inferred globally. Hence, the 
electron densities rie of ~ 4 — 7 x 10"^ cm""^ on small spa- 
tial scales are several times greater than those estimated 
from global radio fluxes, a nd roughly 10 time s those de- 
rived from optical spectra (jlzotov et al.lll999f) . The sen- 
sitive high-resolution measurements presented here are 
able to probe more deeply into the optically-thick star- 
forming region in SBS 0335— 052, and support the idea 
that beam dilution seriously affects low-resolution radio 
observations of such sources. Moreover, optical wave- 
lengths may not be sampling the same regions as these 
radio observations; long wavelengths are better able to 
examine the dusty dense zones in extreme modes of star 



formation such as those in SBS 0335— 052. 

The model fits show that the radio emission appears 
to become optically thick at frequencies ^5-8 GHz. Even 
at 8.5 GHz, the best-fit models suggest that Tff~0.3-0.5. 
Nevertheless, the highest frequency observations pre- 
sented here at 22 GHz are consistent with being optically- 
thin, and we do not believe that the flux densities at 
this frequency have suffered from any significant self- 
absorption. The observed flux densities may, however, 
suffer from losses due to other factors. First, dust within 
the Hii region can absorb ionizing photons, thereby re- 
ducing the amount of ionized gas. Second, if the ISM is 
porous (which we will return to in § 15. ip . a significant 
fraction of the ionizing flux may escape from the region. 

4.3. Ionizing Luminosities 

Thermal radio emission can be used to estimate the 
production rate of ionizing photons from the massive 
star clusters powering these regions, and the ionizing lu- 
minosities can in turn be used to estimate their stellar 
content. 

Following linbS (Il968l ) and lCondonI (|1992f ). 

/ N -0.45 / xO.l 

Qlj,,> 6.32 X 10^2 s-i (0.926) [j^] [^) 

The factor 0.926 corrects for the ~8% helium abun- 
dance in the ionized gas fMez ger fc HendersonI Il967f ) . 
This equation assumes the emission is both thermal and 
optically thin, as both of these criteria are most likely 
to be met at the highest frequency observations, which 
are less likely to contain contaminating non-thermal flux 
and less likely to suffer from self-absorption. We use 
the 1.3 cm flux densities in order to determine Qhyc 
for the thermal radio sources in SBS 0335— 052. We 
have assume d an ionized gas electron temperature of 
re=20000K (jlzotov et al.lll999D . Given that an 07.5V 
star produces an ionizing flux of QLyc= 10"^^ s~^ (here- 
after O*, the equivalent to type 07.5y: ILeithereil il990l : 
lVaccal[l99l iVacca. Garmanv. fc Shul]|ll996n . from Table 
[4] we infer that regions SSCs 1 and 2 contain a minimum 
of approximately 3500 and 2500 equivalent O* stars, re- 
spectively. The entire region S is subject to an ionizing 
flux from the equivalent of on the order ^11,700 O* stars. 

These QLyc values could very well be underestimates 
if a significant fraction of the ionizing flux is absorbed by 
dust or able to escape from the Hii region, possibly due to 
clumping in the ISM, which we believe to be an issue for 
these objects (see § 15.11 and Reines, Johnson., fc HuntI 
l2008f l. Such an effect has been o bserved for Galac- 
tic ultracornpact Hii regions (e.g., iKurtz et al.l Il999l : 
iKim fc Kooll200il ). for which it is estimated that typi- 
cally > 80% of the ionizing flux from the embedded stars 
escapes to an outer diffuse halo. We can estimate an 
upper limit to this effect for sources SSCs 1 and 2 in 
SBS 0335— 052 by assuming that any flux able to leak 
from these regions is still contained in the larger region 
S. This comparison suggests that ~ 50% of the ionizing 
flux may be able to escape the immediate vicinity of SSCs 
1 and 2, and is in reasonable agreement with the escape 
fractio n of 40% estimated bv lReines. Johnson, fc Huri^ 
(|200i) . 



8 



Johnson, Hunt, & Reines 




Fig. 4. — Results from model fits to SSC 1, SSC 2, and Region S (from top to bottom, respectively), (left) Contours of the minima in 
reduced x'^; a clear degeneracy between radius and density is apparent. The diamonds indicate the best-fit purely thermal model, and the 
squares indicate the best-fit model if a non-thermal contribution is allowed, (right) The best-fit model SEDs over plotted on the observed 
data points. Only SSC 1 has a marginally improved fit by allowing for a non-thermal contribution. 



4.4. Age of the SSCs 

The ages of the star clusters m SBS 0335— 052 are not 
weh constrained by our observations. However, the exis- 
tence of thermal radio emission from only SSCs 1 and 2 
suggests that they are younger than the other four SSCs 
detected in the optical with HST. Several other observa- 
tions support t he extreme y o uth o f SSCs 1 and 2. 

For example, IVanzi et"aLl ()2000[ ) observe a very large 
Br7 equivalent width, one of the highest ever observed 
for an extragalactic object. While none of the clus- 
ters show Lya emission, SSCs 4 and 5 are consider- 
ably bright er than SSCs 1 an d 2 in the far-ultraviolet 
continuum (jKunth et al.ll2003l) . This is consistent with 
SSCs 1 and 2 having younger ages and thus suffering 
from more extinction in t he natal cloud. More recently, 
iReines. Johnson, fc HuntI ([2008) find that SSCs 1 and 2 
have ages < 3 Myr based on fitting model SEDs to op- 
tical photometry and measuring Ha equiv alent widths, 
which is consistent with the earlier work of iThuan et al.l 
(|1997f ). We find that the Paa equivalent width is not 



a reliable tracer for age due to the near-IR excess and 
free-free emiss ion that can signifi cantly contaminate the 
continuum (R eines. Johnson, fc H unt 2008). The upper 
limit of ~3 Myr arises because hydrogen emission-line 
widths and stellar colors are insensitive to differences 
in ages less than this, while the ionizing flux is roughly 
constant before the most massive stars have started to 
die. Moreover, the stellar synthesis models are not well- 
calibrated at these extremely young ages. This up- 
per hmit of 3 Myr agrees with previous radio studies 
(e.g. lKobulnickv fc Johnsonlll999[ ) that conclude that the 
SSCs which are detectable as optically-thick thermal ra- 
dio sources are < 1 Myr (based on statistics). SSCs 3-6, 
which are not detec ted in the radio, have ag es ranging 
from ~ 7 — 15 Myr (jReines. Johnson, fc HuntI 12008 1. 

The seemingly insignificant difference in age of a few 
Myr for an SSC could actually be important; ^ 3 — 4 Myr 
is a very "special" epoch in the life of a young star- 
burst. SSCs and their environment undergo a tremen- 
dous amount of development between 1 Myr and 6 Myr, 
making this age range important to target for studying 
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the details of their early evolution. First, this age is just 
after the onset of significant Wolf-Raye t (W-R) popula- 
tions at ~ 3Myr (jLeitherer et al.lll999[) : the W-R com- 
ponent augments by a facto r of ten or so the am plitude 
of the cluster stellar wind ('Leither er et alJ 11992). Sec- 
ond, at ^3.5 Myr, Type II SNe begin to explode, further 
disrupting the ISM through shocks, and increasing still 
more the amount of thermalize d gas streaming out of th e 
confines of the star cluster fe.g. lChevalier fc Clegdll985l ). 
This is the age when the ISM is subjected to an abrupt 
change in mechanical energy input because of the SN 
contribution, and to a lesser extent from the W-R winds. 
The presence of a small, but measurable, W-R popula - 
tion in SSC 3 (iPapaderos et al.ll2006t llzotov et al.ll2006D . 
but not in SSCs 1 and 2, is a confirmation of the younger 
age of these latter clusters. 

During this short period of only a few million years, 
the clusters will transform from being optically obscured 
and extremely luminous in the infrared, to completely 
optically visible with little to no extinction. Very young 
sources are embedded in dense cocoons of gas and dust 
and show rising radio spectra which gradually flatten 
and become fainter as they become increasingly opti- 
cally thin with time ([Cannon k Skillmanl[200l iJohnsoiil 
12004 iHirashita fc HuntI 120061 ). The clusters at ages of 
<4 Myr but older than >3 Myr will develop a wind com- 
ponent in their SEDs, unlike younger clusters without 
W-R stars and SNe. At lower frequencies, the spec- 
tra should begin to steepen around ^S.BMyi because of 
the non-thermal emission from Type II SNe, and finally 
become predominantly non-thermal as the thermal Hii 
emission dies out at ~10Myr. Thus, with sufficient sen- 
sitivity, the radio spectral energy distribution can also 
serve as a diagnostic of a cluster's evolutionary stage. 
As discussed above, the maximum age that SSCs 1 and 
2 could have is ~ 3 Myr, and they could be as young as 
^ 1 Myr. By contrast, SSC s 4 and 5 are ~ 12 - 15 Myr 

eines. Johnson, fc Huntj[2008D . and fully optically vis- 
ible with no detected radio or thermal infrared emission. 

Given t he non-thermal em ission detected on large 
scales by iHunt et"an (j2004l ). it is clear that previ- 
ous episodes of star formation have taken place in 
SBS 0335— 052. However, at the present time we are not 
able to precisely associate the non-thermal emission with 
a specific cluster or clusters. 

4.5. Stellar Content and Star Formation Rates 

The stellar content of the radio Hii regions can be es- 
timated from their Lyman continuum luminosities (ne- 
glecting any leakage or du st absorp tion) using the Star- 
burst99 models of Leither er et al.l ([l999). For a clus- 
ter <3 Myr old, formed instantaneously with a Salpeter 
IMF, 100 Mq upper cutoff, 1 Mq lower cutoff (note that 
reducing the lower mass cutoff to 0.1 Mq increases the 
cluster mass by a factor of ^ 2.5), and 5% solar metallic- 
ity, a IO^Mq cluster has Qlvc^ 7.4 x 10^^ s~^ Assum- 
ing that Qlvc scales directly with the cluster mass, we 
find that the thermal regions around SSCs 1 and 2 are 
powered by stellar clusters with masses ~ 5 x 10^ Mq 
(Table 3]). Since the ionizing luminosity is roughly con- 
stant from 1-3 Myr, this result is independent of age, as 
long as the age of the clusters is <3 Myr. These masses 
are roughly a factor of two low er than those found by 
iReines. Johnson, fc HuntI (|2008f) . because of the under- 



luminosity of the radio emission relative to the optical 
SEDs (see gSI]). 

In these young massive clusters, the number den- 
sity of just the O* stars is high: ~3000 such stars in 
a sphere of 4pc radius gives ~llpc~^. These densi- 
ties are such that the extent of a typical stellar wind 
would significantly overl ap with that of it s neighbor 
in less than IMyr (e.g., IWeaver et"all I1977D . The ef- 
fect of th ese "colliding winds" should be strong X-ray 
emission ()Pittard fc Stevens! [l997t iMvasnikov fc Zhekovl 
,1993) . over and above the diffuse X-ray emission ex- 
pected from the adiabatic interaction of the stellar wind 
(jChevalier fc Clegd [19851 iCanto et al.l[2000l ). Compact 
X-ray emission is indeed observed in SBS 0335— 052 
([Thuan et al.l [200 ^ , although the spatial resolution of 
those observations is insufficient to associate it only with 
SSCs 1 and 2. 

At an age of <3Myr, the individual clusters SSCs 1 
and 2 would together provide a bolometric luminosity 
of Lboi~ 1-5 — 1.9 X 10^ Lq, according to the mod- 
els of Starburst99 for an instantaneous burst with a 
metallicity of 5% solar. This is in reasonable agreement 
with the infrared luminosity of SBS 0335—05 2 estimated 
from the dust emission, 1.4 — 1 .5 x 10^ Lq ([Hunt et al.l 
I2005al: lEngelbracht et all [20051 ). The additional ioniz- 
ing luminosity in the entirety of region S would yield 
-^boi'^ 3 x 10^ Lq. This excess at larger scales could 
mean that the dust is intercepting only about 50% of the 
ionizing radiation globally. Consequently, the bolometric 
luminosity inferred from infrared observations could be 
underestimated. 

The instantaneous star formation rate (SFR ) can be es- 
timate d using the QLyc luminosity. Following iKennicuttI 
([199^) : 

SFR(Afo year^i) = 1.08 x lO"'"^^ QLyc{s~^)- (7) 

The inferred QLyc value from the 1.3 cm emission yields 
a SFR~ 1.3Moyr~^ for region S (270 pc in diameter), 
which is quite substantial. Alternatively, with an upper 
limit to the age (<3 Myr) and an estimate of the stellar 
mass of the SSCs 2x5x10^ Mq, corrected by a factor 
of 2.5 to extend the IMF down to O.IM0 and dividing 
by a factor of 0.6 to roughly correct for the fraction of 
escaping radiation), we can simply determine an aver- 
age SFR, and find ^1.4 M , Tiyr~^ . By comparison, using 
Ha emission. iThuan et al.l ( 19971 ) find a SFR for the en- 
tire galaxy of ^ 0.4Mq yr~ ^ , already among the highest 
SFRs observed in BCDs ([Fanelli. O'Connell. fc Thuad 
|1988() . This discrepancy may suggest that a substantial 
portion of the current star formation that is traced by 
thermal radio emission in SBS 0335—052 i s hidde n from 
view in optical recombination lines (H unt et al.l 12001). 
However, there is at least a ~ 30% dispersio n between 
SFR calibrations from different authors (e.g. IKennicuttI 
Il998l and references therein), which could partially ac- 
count for the observed discrepancy. This is similar to 
the difference in star formation derived fr om optical and 
radio observations for the BCD Haro 3 ([Johnson et al.l 
|2004[ ). Clearly the nascent star formation in region S is 
a major event for this dwarf galaxy, and serves to high- 
light the impact a small region of star formation can have 
on a small dwarf galaxy. 

The relevance of such an impact can be perhaps bet- 
ter appreciated by considering that this relatively high 
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(for a dwarf galaxy) star formation rate takes place in 
a very small region. Considering a volume with radius 
R ~135pc (as for Region S), we would derive a star- 
formation per unit area of ~23 M0yr~^kpc~^. This is 
comparable to, although lower than, the st arburst inten- 
sity li mit of 45 Mq yr^^ kpc~^ derived bv lMeurer et al.l 
(|1997f ) from observations of starbursts at various wave- 
lengths. Star formation in the SSCs in SBS 0335—052 is 
occurring very close to the maximum intensity observed 
in the universe, and, unlike larger systems, this mode of 
star formation dominates the energy output of the entire 
galaxy. 

5. THE STARBURST IN SBS 0335-052 

In what follows, we discuss the consequences of our re- 
sults in the context of the galaxy itself, and their impact 
on our knowledge of how stars form in low-metallicity 
environments. 

5.1. The Porous Circum- cluster Medium 

Several different results suggest that the ISM surround- 
ing the natal clusters in SBS 0335— 052 is porous and 
clumpy. First, the spectral fit of Region S suggests that 
there is ionizing radiation beyond the strict boundaries 
of the SSCs 1 and 2. The volume emission measure V 
and the number of ionizing photons QLyc (see ^4.3p of 
region S are nearly twice as large as the sum of SSCs 1 
and 2 (Tabled]). This implies that there are ionized zones 
outside of the strict confines of the SSCs that contribute 
to the radio emission. The extended emission could be 
caused by a low level of distributed star formation. Al- 
ternatively, it could be caused by ionizing radiation that 
has escaped from the compact Hii regions via a porous 
ISM. This effect due to ionizing radiation leakage has 
also been o bserved for u ltra co mpact Hil region s in the 
Milky Way ()Kurtz et al. 1999: Kim fc Kool[200l[ ). 

In previous work, iThuan et al.l ^2005) inferred 
a clumpy ISM in SBS 0335-052 based on their 
H2 absorption-line obs e rvatio ns. More recently, 
iReines. Johnson, fc HuntI (|2008D also reveal evidence for 
a porous and clumpy ISM surrounding the young SSCs 
1 and 2: the measured ionizing luminosities from Ha, 
Paa, and optically-thin free-free radio emission are lower 
than expected compared to the optical SEDs. Reines et 
al. infer that ^ 40% of ionizing photons from the stellar 
continuum are leaking out of the immediate vicinity of 
the clusters before intercepting hydrogen to ionize. Thus, 
a significant fraction of the ionizing flux is escaping with- 
out contributing to the measured ionized gas emission. 

Finally, despite the substantial extinction de- 
riv ed fro m mid-IR observati ons toward this region 
(iThuan et all 11 999: Pla nte fc Sauvagd 120021 : iHunt et all 
I2005at iHouck et al.. ,200ir some UV con tinuum light 
is still able to escape (jKunth et "all 120031 ) . Moreover, 
optically-derived extinctions are quite low (|Izotov et al.l 
[1997; Reines, Johnson, & Hunt 2008). This would natu- 
rally follow from inhomogeneities in the surrounding dust 
that would simultaneously allow for large extinctions de- 
rived from mid-IR emission originating from dense dust 
clumps (as opposed to a uniform screen) and low extinc- 
tion regions from re latively transparent lines-of-sigh t into 
the birth cloud f see IReines. Johnson, fc Huntll2008l for a 
full analysis and discussion). 

Given the impact that clumping of the ISM can have 



on the resulting observed spectral energy distribution, as 
seen in this case for SBS 0335—052 it is important to keep 
the potential effect of dumpiness and inhomogeneities in 
mind when interpreting similar data. 

5.2. Toward Reconciling the Inferred Extinction 

The nature of the observed and predicted line ratios in 
SBS 0335 —052 has b een a source of discussion in the liter- 
ature fe.g.lHunt et a l. 2001, 2004; Thompson et al. 200^; 
IReines. Johnson, fc HuntI 12008: Tho mpson et al.l 120081. 
In particular, an apparent discrepancy between the ther- 
mal radio emission and infra red line fluxes was first dis- 
cussed by iHunt et al.l (I2004D using global fluxes, and by 
iThompson et al.l (|2006f using a preliminary version of 
the radio obser vations no w published here. The intrinsic 
thermal flux in iHunt et al. (2004) , inferred from fitting 
the radio SED, is about 50% lower than we find with 
these new observations, and thus appeared to be under- 
luminous with respect to previous near-infrared (NIR) re- 
combination line measurements (e.g.. B ra and Br7) . The 
unpub lished radio flux densities used bv lThompson et aTl 
(|2006f) were derived from a smaller aperture than used 
in this paper, consequently resulting in lower flux densi- 
ties. Using the updated values, we find that the observed 
and predicted Paa values (1.77 x 10"^'* and 4.05 x 10"^"* 
erg s^^cm"^, respectively) are in good agreement for a 
native extinction of Apaa = 0.9, which is significantly 
lower th an the Apaq = 1.64 adopted b y Thompson et 
al. (see IReines. Johnson, fc Hunti l2008l for further dis- 
cussion). Thus, we conclude that observed emission line 
fiuxes and line ratios are fully consistent with the fiux 
densities obtained from the new radio observations pre- 
sented here. 

5.3. Stellar Winds? 



iThuan fc Izotovl (|1997( ) found evidence for P Cygni 
profiles in SBS 0335-052, integrated over the 2x2 
arcsec^ aperture of the Goddard High-Resolution Spec- 
trograph aboard the HST This aperture is sufficiently 
large to encompass all 6 SSCs in SBS 0335—052, so it was 
difficult to pinpoint their origin. The most likely scenario 
seemed to be that they arise from mas sive stellar winds 
in th e older clusters to the northwest (|Thuan fc Izotovl 
Il997t l. rather than the younger ones to the southeast. 
This conclusion is supported by later observations of 
wider Ha pro files i n SSCs 4 and 5 than in SSCs 1 and 2 
(|lzotov et al.ll2006l). 

More recentlv. IThompson et al.l (|2006f ) interpreted the 
radio fiux densities of SSCs 1 and 2 as resulting from 
a dense stellar wind causing self-absorption at all ob- 
served radio wavelengths. If there were such a wind 
powering the radio emission in SBS 0335—052, we would 
first expect no spectral turnover at the highest fre- 
quencies observed, in contrast with the spectral energy 
distribution of region S. We would also expect signifi- 
cant NIR recombination line excesses relative to "nor- 
mal" Hii-region emission. Infrared recombination lines 
are expected to show the effect of dense stellar winds 
much more than the radio, because the radio signa- 
ture of such winds is at least two orders of magnitude 
weaker t han the effect on NIR recombination line opac- 
ity (e.g. ISmith et al.l 119871 : ISimon et al.l ll983l Such an 
effect is observed for the massive O stars in t he Galactic 
Center Arches Cluster (|Nagata et al.l 119951 : iLang et aLl 
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l2001f) . Relative to the ratio in Hii regions, the IR- 
line-to-radio-continuum ratio for those stars is >250; in 
other words, the IR line emission generated from ion- 
ized winds around massive stars is roughly 250 times 
stronger relative to the radio continuum than it would 
be in an Hii region. This is entirely consistent with stel- 
lar wind theory (e.g..lKrohk fc Smithlll981l : ISimon et all 
flOSllSmith et alJIlOSTl) . However, NIR line excesses are 
not observed in SBS 0335— 052^. Because of the consis- 
tency with our new ra dio observations and the NIR re - 
combination lines (e.g.. iReines. Johnson, fc Huntll2008D . 
we conclude that SSCs 1 and 2 do not show evidence 
for massive stellar winds. The lack of observed winds 
is consistent with both their young age of <3Myr (no 
W-R population) a nd low mctallici ty (mass loss goes 
roughly as seej Abbott. 1982; Leitherer et al.ll 19921 : 

ICrowther et al.|[2002[ ). Because the combined effects of 
young age and low metallicity conspire to reduce wind 
power, a lack of winds alone does not constrain the ages 
of SSCs 1 and 2 to be young. However, winds have been 
observed in other clusters in SBS 0335— 052 (see above), 
which suggests that the difference could be the older age 
of those SSCs. 

5.4. Dearth of Non- Thermal Emission at High 
Resolution 

The data presented here are fully consistent with 
purely thermal models. The results of non-thermal mod- 
els only give a marginally better fit for source SSC 1 
than the purely thermal models, and at the cost of ad- 
ditional free parameters. Nevertheless, if non-thermal 
emission is included in model fits, SI is consistent with 
having a non-thermal fraction at 5 GHz of ^ 5% with a 
non-thermal slope of —0.8. The best fits for SSC 2 and 
Region S are purely thermal. In order to include a suffi- 
ciently weak non-thermal (synchrotron) emission in the 
model fits, these data could at most support roughly a 
single SN observed soon after explosion, or alternatively 
a larger number of SNe at somewhat later epochs. 

None of the SSCs in SBS 0335-052 have detectable 
non-thermal radio emission in the data presented in this 
paper. At best, SSCs 5 and 6 have marginal radio de- 
tections at 1.3 cm. SSCs 4, 5, and 6 all appear to be 
sufficiently old that they are likely to contain a number 
of SNR. H owever, based on a comp arison with the pre- 
dictions of iHirashita fc Hun^ (|2006D . the expected non- 
thermal flux originating from these clusters is slightly 
below the sensitivity limits of these data, which are rela- 
tively high-frequency and not optimized for the detection 
of non-thermal sources. 

//the observations were to definitively detect the pres- 
ence of such non-thermal emission from these natal clus- 
ters (assuming it were not due to a projection or confu- 
sion effects), it would serve as a strong lower limit on the 
age of the star cluster. Stellar evolutionary models pre- 
dict that no SNe cou ld have exploded before ~3.5Myr 
(|Leitherer et al.|[T999f ). This age would be in co n tradic - 
tion to the results from Kobul nickv fc JohnsonI (jl999) , 
which suggest (based on statistics) that the lifetime of 
the "ultradense Hii regions" is < 1 Myr. In addition, 

6 IHunt et al] 1I2OOI ') found such an excess, but their thermal 
radio flux is half of what we deduce from the new high-resolution 
observations presented here. 
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Fig. 5. — A comparison between the flux densities observed 
at high-reso lution (this paper) and the flux densities observed by 
IHunt et al] g004), which are sensitive to larger spatial scales. The 
different symbols (open circle, star, square, and triangle) corre- 
spond to the global (low-resolution) data, Region S, SSC 1, and 
SSC 2, respectively. The best-fit models are also shown: the shape 
of the global spectrum is indicative of a strong non-thermal com- 
ponent. 

upper fimits o n the ages of SSCs 1 and 2 of < 3 Myr were 
determined bv lReinesT Johnson, fc HuntI (|2008f ) usmg op- 
tical and infrared observations; these results provide ad- 
ditional support for the young ages of these clusters, and 
thus the expected lack of non-thermal emission. 

These results may app ear to be i n co ntradiction 
with previous results from iHunt et~al\ (|2004l ) , in which 
a non-thermal c omponent was detected. However, 
iHunt et"al] (|2004[ ) obtained lower resolution observations 
of SBS 0335—052 that are sensitive to larger-scale spatial 
structure than the observations presented here (by virtue 
of interferometers acting as spatial filters) . The flux den- 
sities derived from the lower resolution observations are 
over-plotted on the high-resolution observations in Fig- 
ure [51 The Hunt et al. data clearly indicate the presence 
of non-thermal emission in SBS 0335— 052, as manifest 
in the negative slope of the spectral energy distribution. 
A comparison between the two data sets demonstrates 
that the observations presented here are resolving out 
a substantial fraction of the longer waveleng th emission 
at 6 cm (roughly 44% of that observed by iHunt et al.l 
[200l . while recovering all of the emission at shorter 
wavelengths. The unresolved 20 cm source in the lower 
resolution observations is also offset to the northwest of 
SSCs 1 and 2 (see Figure [3]), but there is no evidence in 
those observations for a peak in non-thermal emission at 
the location of the thermal radio sources. However, the 
center of the diffuse non-thermal emission could possibly 
coincide with a faint optical source, visible in the ACS 
image toward SSC 3. 

We ascribe this difference to the diffuse nature of 
synchrotron halos around massive st ar clusters (e.g., 
'Cann on et all 120051 : iHunt et al.|[2005bD . The dearth of 
non-thermal emission in these high-resolution observa- 
tions suggests that the non-thermal emission is diffuse 
in nature, and not obviously associated with the com- 
pact thermal radio sources we detect. If the synchrotron 
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halo of IZwl8 were to be placed at the distance of 
SBS 0335— 052, it would be resolved out by our high- 
frequency observations (LAS~2"). Such halos are domi- 
nant at low frequencies and have low surface brightness, 
which makes them difficult to detect with the observa- 
tions presented here. Hence, we conclude that our obser- 
vations are not inconsistent with previous ones. More- 
over, the presence of non-thermal emission in close vicin- 
ity to the compact thermal radio sources also suggests 
that significant star formation has taken place preceding 
the birth of sources SSCs 1 and 2. 

5.5. The Importance of Metallicity in Massive Star 
Cluster Evolution 

Natal super star clusters presumably emerge from their 
birth cocoons via a combination of expansion of the high 
pressure regions, radiative dissociation of molecules and 
dust, and stellar winds. The interplay of these processes 
will affect the timescales in the evolution of natal clus- 
ters, and in particular how rapidly a cluster emerges from 
its birth material. Each of these processes has a de- 
pendence on metallicity, and thus the low metallicity of 
SBS 0335—052 could potentially have a strong combined 
effect on the early evolution of the natal clusters that it 
hosts. 

First, the pressure of the ionized gas within the cluster 
may be affected by metallicity. Metallicity has a fun- 
damental role in the rapidity of radi ative cooling of the 
ISM within and outside the cluster (iMcCray &: KafatosI 
[19871 : iSihch eTalllM)! pTeno rio-Tag le et al.ll2005ar ). and 
therefore the pressure of the region will be affected. 
The average pressure of the ionized gas in the clus- 
ters implied by the densities we infer from the spec- 
tral fits is quite high, ^ 1.4 x 10~* dynecm"^. Given 
the ionized gas electron temperature of Te=20000K 
(jlzotov et a l. 1999) is roughly twice the "typical" value 
expected for Hii regions, it follows that the pressure 
would be correspondingly twice as high as that found 
in a typical Hii region. A comparison with other galax- 
ies hosting natal clusters suggests that the inferred pres- 
sure for the SSCs 1 and 2 in SBS 0335-052 is simi- 
lar to in other SSCs in dwarf galaxies (e.g., He 2— 10, 
IC4662, NGC5253, IIZw40 iJohnson fc Kobulnick^ 

et al .l l200flHBeck et al 



2003H Johnson et al.ll2003H TiI r 



mm . However, given the uncertainties in the model fits 



to the spectral energy distributions, the inferred pres- 
sures are, at best, accurate to within a factor of two. 
More precise measurements of the pressure in natal SSCs 
from a range of environments will help to address this is- 
sue. 

The metallicity can also affect the abundance of 
molecules and dust in the natal cocoon. First, at 
low metallicities, the relative amount of molecules and 
dust is expected to be lower than in higher metallic- 
ity counterparts. Moreover, the dust itself has a criti- 
cal role in the formation and survival of molecules (e.g . 



HoUenbach et al.''l971; Salp eteij[T977l: ISavage et anil977t 
HoUenbach fc Tielens .1997,) , thus there is an intimate 



connection between metallicity, dust, and molecules. 
In addition, lower metallicity systems will also have 
harder radiation fields and create a more hostile environ- 
ment for molecules and dust (e.g. iMadden et al.l l2006t 
iLebouteiller et al] 120071 ). These effects would tend to 
shorten the timescale for a natal SSC to emerge from its 



birth material. SSCs 1 and 2 in SBS 0335-052 are al- 
ready associated with some optical emission, consistent 
with this scenario, and it is possible that the low metallic- 
ity environment has accelerated the emergence process. 
However, the partial emergence of these clusters could 
also simply reflect their nominal evolutionary state, and 
we are not able to discriminate between these possibili- 
ties with existing data. 

Another mechanism by which metallicity may affect 
the early evolution of massive clusters is via their stel- 
lar winds. It is well known that the eff ect of massive 
stellar winds decreases with metallic i tv (lAbbottl IT981: 
Leitherer et al.lll992Hyink et al.l[200lHKudritzki fc Puld 



2000t iCrowther et al.l |2002| ). primarilv due to reduced 
line blanketing. Even before the onset of the W-R phase 
and the first supernova in a cluster, the mechanical en- 
ergy input from the massive star winds alone is more than 
an order of magnitude weaker for systems with metallic- 
ity as low as S BS 0335—052 compare d to solar metallicity 
counterparts (jLeitherer et al.l[l999l ). In contrast to the 
effects discussed above, the reduced mechanical energy 
input for a low metallicity cluster would tend to lengthen 
the embedded phase - radiation pressure would not be 
able to clear away the surrounding material as effectively. 
However, at the present time, it is unclear which (if any) 
of these effects will have a dominant role in the early 
evolution of massive star clusters with low metallicity. 

6. SUMMARY 

This paper examines the radio continuum properties 
of current star-forming regions in the extremely low- 
metallicity galaxy SBS 0335—052 using high-resolution 
VLA observations. Two intense star-forming regions are 
detected as luminous thermal radio sources that appear 
to be extremely young SSCs. The total ionizing flux 
of the southern star-forming region in the galaxy sug- 
gests the presence of ~ 12, 000 "equivalent" 0-type stars, 
and the inferred instantaneous star formation rate for the 
radio-detected natal star clusters alone is ~ 1.3 Moyr"^ 
or ~ 23 Moyr^^kpc"^, which is close to the starburst 
intensity limit of 45 MQyr~^kpc~^. This star formation 
rate derived from thermal radio emission also suggests 
that previous optical recombination line studies are not 
detecting a significant fraction of the current star forma- 
tion in SBS 0335-052. The observations presented here 
also suggest that up to ^ 50% of the ionizing flux could 
be leaking out of the compact Hii regions; this in is agree- 
ment with previous work that suggests the interstellar 
medium surrounding the natal clusters in SBS 0335—052 
is porous and clumpy. Model flts to the radio spectral 
energy distribution indicate that the global mean density 
in the youngest SSCs is Ug > 10^ — 10** cm~^. 

A compa ri son b etween these radio data and those of 
iHunt et all (|2004[ ) suggest that the observations pre- 
sented here have resolved out a significant contribu- 
tion from non-thermal emission over spatial scales not 
detected in these high-resolution observations. The 
presence of this diffuse non-thermal emission may in- 
dicate a synchrotron halo in the southern region of 
SBS 0335— 052, likely associated with a previous episode 
of star formation in proximity to the natal clusters de- 
tected here that has not yet been identified. Such a pre- 
vious generation of massive stars will have a rapid impact 
on the chemical enrichment in this very low-metallicity 
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galaxy. 

A fundamental question is how SSCs form and evolve 
at low metallicities. While we cannot presume to an- 
swer this with these observations alone, it is clear that 
SBS 0335— 052 is undergoing an extreme burst of star 
formation, the results of which are likely to significantly 
impact the galaxy's future evolution. The low metallie- 
ity in this galaxy may affect the evolution of the natal 
clusters in a variety of ways, including inefficient cool- 
ing, hard radiation fields, and weak radiation pressure 
from line blanketing. However, a larger sample of low- 
metallicity natal clusters with panchromatic observations 
will be required to address these issues. These observa- 
tions do, however, conclusively demonstrate that massive 



star clusters are indeed able to form in low-metallicity 
environments, which is an important step toward our 
understanding of globular cluster formation in the early 
universe. 
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